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Abstract 

We present a preliminary measurement of the branching fraction for the decay of the neutral 
B meson into the final state J/07r + 7r~. The data set contains approximately 56 million BB pairs 
produced at the T(4S) resonance and recorded by the BABAR detector at the PEP-II e + e~ collider. 
The result of this analysis is B(B° —* J/ipir + ir^) = (5.0 ± 0.7 ± 0.6) xlO -5 , where the first error is 
statistical and the second is systematic. 
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1 Introduction 



In this paper, we report a preliminary measurement of the B° — > J/ip Tr + ir~ branching fraction. The 
motivation for studying these decays comes in part from the possibility of including the J/i/j p° mode 
in the BABAR sin2/3 analysis Q. It is also expected that due to the Cabibbo and colour suppression 
of these decays, interference with rare or exotic processes, such as a box diagram containing charged 
Higgs bosons, could be significant. These effects may appear as deviations of the branching ratio 
from the Standard Model prediction of B(B° -> J/^ir + ir~) = (4.6 ± 0.8) x 1(T 5 [§. The only 
previous measurement related to this analysis is an upper limit on the B° — ► J/ipp° branching 
fraction by CLEO §, B(B° -» J/^p°) < 2.5 x 10" 4 at 90% confidence level. 

The method used is a fit to the distribution of the invariant mass of the two pions, M(ir + Tr~), 
with the goal of isolating the separate components. Both the B° — > J/tp p° and B° — > J/ipir + ir~ 
(non-resonant) contributions are modelled in the fit function. However, in the present analysis the 
branching fraction is quoted for the sum of these two modes, as the individual fractions are highly 
correlated and depend sensitively on model assumptions. The total yield is only sensitive to the 
normalization of various background sources and so is quite stable. It is expected that separation 
of the two components will be possible with a considerably larger data set, when a mass dependent 
angular analysis of the decay products will be possible. 

2 The BABAR detector and data set 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. 
The BABAR detector is described in detail elsewhere ||. Surrounding the beam pipe is a silicon 
vertex tracker (SVT) to provide precise measurements of positions and angles of charged particles 
emerging from the interaction region and stand-alone track finding, particularly for particles with 
momentum below 120 MeV/c. Outside this is a 40-layer drift chamber (DCH) filled with an 80:20 
helium-isobutane gas mixture to minimize multiple scattering. The DCH provides measurements of 
track momenta, as well as energy-loss measurements that contribute to charged particle identifica- 
tion. Surrounding the drift chamber is a novel detector of internally reflected Cherenkov radiation 
(DIRC) that allows charged hadron identification in the barrel region. Outside the DIRC is a 
CsI(Tl) electromagnetic calorimeter (EMC), which is used to detect photons and neutral hadrons 
and to provide electron identification. The EMC is surrounded by a superconducting coil, which 
provides a magnetic field of 1.5T for momentum measurements. Outside the coil the flux return 
is instrumented with resistive plate chambers interspersed with iron (IFR) for the identification of 
muons and long-lived neutral hadrons. The GEANT4 Q software is used to simulate interactions 
of particles traversing the BABAR detector. A coordinate system is defined with the z axis aligned 
along that of the detector in the electron beam direction. 

The data sample used for the analysis contains approximately 56 million BB pairs, corre- 
sponding to 51.7 fb _1 taken near the T(4S) resonance. An additional 6.4 fb -1 of data, taken 
approximately 40MeV below the T(4S) peak, was used in studies of the effect of light qq pair 
backgrounds. 

3 Analysis method 
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3.1 Event selection and B reconstruction 

The aim of the event selection is primarily to accept events containing a B° — > J/ipir + ir~ decay 
while rejecting background both from u, d, s, c quark continuum and T(4S) events which do not 
contain a B° — > J/tpir+ir^ decay. 

Events containing BB pairs are selected based on track multiplicity and event topology. Tracks 
in the polar angle region 0.41 < 9^ < 2.54 rad, with transverse momentum greater than 100 MeV/ c 
are required to pass quality cuts, including number of DCH hits used in the track fit and impact 
parameters with respect to the nominal beam spot in the r — <f> and r — z planes. At least three 
tracks must pass the above selection. To reduce continuum background the ratio of second to zeroth 
Fox- Wolfram moments ||, R2 = H^/Hq, is required to be less than 0.5. The total charged and 
neutral energy must be greater than 4.5 GeV in the fiducial region of the detector (charged tracks 
in the DCH and neutral clusters in the EMC) and the primary vertex of the event must be within 
0.5 cm of the average measured beam spot position in the plane transverse to the beamline. 

Electron candidates must satisfy the requirement that the ratio of calorimeter energy to track 
momentum lies in the range 0.75 < E/p < 1.3, have a cluster shape and size consistent with an 
electromagnetic shower, and DCH dE/dx consistent with an electron. If an EMC cluster close to 
the electron track is consistent with originating from a bremsstrahlung photon, it is added to the 
electron candidate. 

Muon candidates must satisfy requirements on the number of interaction lengths of IFR iron 
penetrated of N\ > 2, the difference in the measured and expected interaction lengths penetrated 
of \N\ — N^ xp \ < 2, the position match between the extrapolated DCH track and the IFR hits, and 
the average and spread of the number of IFR strips hit per layer. 

Pion candidates are accepted if they originate from close to the beam spot and fail a charged 
kaon identification algorithm that is designed to reject pions. The algorithm uses dE/dx information 
from the SVT and DCH, and Cherenkov angle and number of photons from the DIRC. 

Tracks are required to lie in fiducial ranges within which the efficiency is well known from control 
samples and the material in the detector is accurately modelled in the simulation. The accepted 
ranges in laboratory polar angle are 0.41 < O^b < 2.41 rad for electrons, 0.3 < #i a b < 2.7 rad 
for muons, and 0.35 < #i a b < 2.5 rad for pions. These ranges correspond approximately to the 
geometrical acceptances of the EMC, IFR and DIRC, respectively. 

J/ip — > e + e _ or /^ + /U~ candidates are formed from pairs of identified electrons or muons which 
are fitted to a common vertex and must lie in the invariant mass interval 2.95 (3.06) to 3.14 GeV/ c 2 
for the e + e~ (// + //~) channel. 

B° candidates are formed by combining a J/tp candidate with pairs of oppositely-charged pion 
candidates that are fitted to a common vertex. We require that both two-prong vertices coincide 
spatially by demanding 

> ^ — — <- (\ (\\ 

'(^) 2 +«) ' [) 

where x l and <r(x l ) are the positions and errors of the vertices in the ith spatial dimension. Further 
selection requirements for B° candidates are made using two kinematical variables: the differ- 
ence AE between the energy of the candidate and the beam energy E^^ in the center-of-mass 



frame and the beam-energy substituted mass m^s = y (-E'beam) 2 — (Pb^) 2 - After applying the loose 
requirements 5.2 < ttt-es < 5.3 GeV/ c 2 and \AE\ < 120 MeV, if more than one B° candidate re- 
mains in the same event only the candidate with the smallest \AE\ is kept. Those candidates 
with |mEs — 5279. 0| < 9.9 MeV/ c 2 and \AE\ < 39 MeV form the final signal sample. These ranges 
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correspond to 3a in the expected resolutions for m^s an d AE. After all selection criteria have been 
applied, there are 212 events remaining. 



3.2 Fitting to the tt + ti mass spectrum 

An unbinned, extended maximum likelihood fit is performed on the invariant mass distribution of 
the two pions for selected events. We consider five categories of events: 

• B° — » J/i/jtt + it~ (non-resonant) events 

• B° — > J/ip p° events 

• B° -> J/ipK° s {K° s -» vr+vr") events 

• Background from events with a fake J/ip 

• Inclusive- J/ip background (that originating from events containing a real J/ip ) 

A probability density function (PDF) is constructed for each of these five cases. The total PDF is 
then formed from the sum of the five PDFs and fitted to the data. The B° — » J/ip K® mode is not 
considered signal for the purposes of determining the branching fraction B° — > J/ipir + ir~. 

3.2.1 The signal and K Q S PDFs 

The PDF used to model the B° — > J/ip 7r + 7r~ (non-resonant) mode is the two-body phase space 
distribution of the pion pair and is given by 

Fn.rXm) = {m- M on ) A x (M off - m) B , (2) 

where the kinematic limits are M on = 2M m and M Q fj = M B o — Mji^. The parameters A,B are 
determined by fitting to the M(7r + 7r~) distribution in simulated B° — > J/ip ir + ir~ (non-resonant) 
events. 

The PDF used to model the B° — > J/tp p° mode is a relativistic p-wave Breit-Wigner @ function 
multiplied by the phase-space shape described above. The p° mass and width for the Breit-Wigner 
are fixed to their PDG values ||. 

The B° — > J/ip Kg two-pion invariant mass distribution is modelled by a single Gaussian function 
with the mass and width fixed to values obtained by fitting a sample of simulated J/ip Kg events. 



3.2.2 The background PDFs 

The M(tt + tt~) distribution for events without a real J/ip is derived from a fake- J/ip sample from 



data. The sample is selected as in Sec. |3.1| except that at least one of the lepton candidates must 
fail the appropriate particle identification requirements, thus providing a high statistics data set 
that models the shape of the non-J/ip background. In addition the final kinematic requirements 
are relaxed to 5.2 < rngs < 5.3GeV/c 2 and \ A.E\ < 120 MeV, which further enlarges the sample. A 
Monte Carlo study confirms that the M(tt + tt~) distribution obtained with this procedure correctly 
describes the shape of the non-J/0 background. The resulting distribution is parameterized using 
the sum of two Weibull functions (see Eq. ^) and a Breit-Wigner. 

The M(tt + 7t~) background shape from events containing a real J/ip is obtained with a sample 
of simulated B — > J/ipX events equivalent to a luminosity of 81 fb _1 . Events where the system X 
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is 7r + 7r (non-resonant), p° or K® (ir^ir ) are removed from the sample. The resulting shape is 
well described by a Weibull function ||: 

%H = CV(m- M on )( c -^ x exp[-U(m - M max f] 

V = (C-l)/(C(M max -M on f), (3) 

where M max (the peak position) and C are determined from the fit to the simulated events and 
M on is the kinematic limit as described above. 



3.2.3 Background normalization 

The normalization of the background components is obtained from control samples in data. The 
level of non- J/ip background is obtained from the data J/ip -sideband sample. This sample is selected 
by accepting only events where the J/ip — > e + e~ candidate lies in the invariant mass region 3.156 
to 3.300 GeV/c 2 or the J/ip -> fj, + fi~ candidate lies in the region 2.980 to 3.024 GeV/c 2 or 3.156 to 
3.300 GeV/c 2 . The rriES distribution is fit to an ARGUS function Q for B candidates in which the 
J /if) — > e + e~ or J/ip — > n + fi~ candidate lies in the J/ip sideband. The fit is done separately for 
e + e~ and /J- + n~ candidates to obtain the number of background events passing the final kinematic 
selection in each case. These numbers are then reweighted to make them correspond to the fitted 
level of background in the M{e + e~) or M(/x + n~) distribution in the J/ip mass interval. 

The normalization of inclusive- J/ip background is obtained from the distribution of rngg for 
events in the AE signal region. The rriES distribution is parameterized by a Gaussian (to represent 
signal and peaking background) and an ARGUS function, which has one shape parameter and an 
endpoint fixed by the average beam energy in the center-of-mass frame. The peaking background 
is defined as that which accumulates at tues = 5.279 GeV/c 2 . It originates from B — ► J/ipX decays 
such as B — > J/ipK* . The non-peaking component of the inclusive- J/ip background comes from 
subtracting the non-J/ip contribution, on the basis of the scaled sideband events described above, 
from the total ARGUS background in data. We then add a peaking component, which is obtained 
from the Gaussian part of the rriES distribution in B — > J/ipX simulation (with X = ir + ir~ (non- 
resonant), p° and K®(ir + ir~) removed). The total background from events containing a real J/ip is 
then the sum of peaking and non-peaking components. The peaking component comprises only 5% 
of the total background. Thus any associated uncertainties, such as those on the branching ratios 
used in the J/ipX simulation, will have a relatively small effect on the final systematic error. 



3.3 Branching fraction determination 

The branching fraction is obtained from 



N B o x e nn x B(J/ip — ► £+£-) 

where A^^,,- is the yield of resonant and non-resonant components obtained from the fit, iV B o is 
the total number of B° and B° in the data sample [|l(]] and e W7V is the signal efficiency. The J/ip 
branching fraction B(J/ip — > £ + £~) is fixed to the current world average value @. We assume that 
the branching fraction T(4S) -» B°B° is 50%. 

The signal efficiency in Eq. |I| is derived from simulated signal events. Imperfect modelling of 
particle identification is corrected by studies of independent control samples derived from data. 
The lepton and pion particle identification efficiencies are measured in data with control samples of 
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p + p 7, p + p e+e , e + e , e + e 7, D* + — > D°7r + (D° — > i^T 7r + ) and — > 7r + 7r . The efficiencies 
are determined as a function of momentum, and polar and azimuthal angle. We find e(J/ipp°) = 
(27.0 ± 0.3)%, e(J/^vr + vr-, non-resonant) = (26.5 ± 0.3)% and e(J/ipK° s ) = (2.31 ± 0.06)%. The 
final corrected signal efficiency is taken as the average of the Jfyp and J/07t + 7t~ (non-resonant) 
efficiencies and is found to be (26.8 ± 0.2)%, where the error is from Monte Carlo statistics. 



4 Results 



The full model for the M(ir + ir~) mass distribution is obtained by summing the five PDFs described 
in sections 3.2.1 and |3.2.2| . A likelihood fit is performed on the M(tt + tt~) distribution in data with 
the normalization of the non-J/0 background fixed to 36.9 events and the inclusive- J/0 background 
to 55.1. Thus the only parameters that are allowed to float in the fit are the numbers of J/07r + 7r~ 
(non-resonant), J/ip p° and J/ipK® events. The result of the fit is shown in Fig. [l] where the data 
has been binned and overlayed. The x 2 °f the curve and the data as binned in Fig. [l] is 40.9 for 
38 data points. The fit finds N wn = 89.9 ± 13.3 events, where 43.3 ± 12.7 are included in the p° 
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Figure 1: Distribution of the invariant mass M(7r + 7r~) for events passing all selection criteria. 
The solid line is the result of the unbinned likelihood fit. The dotted (dashed) line represents the 
background from non-J/ip (inclusive- J/0 ) events. 



component and 46.6 ± 14.5 are included in the non-resonant component. The number of events in 
the Kg component is 29.0 ± 5.5. Inserting the result into Eq. [| yields the branching ratio 

B(B° -» J/07T+O = (5.0 ± 0.7) x HT 5 , 
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Table 1: Number of events for each signal and background category floated in the fit. 
J/ifj7r + ir~ (counting) is the J/07r + 7r~ yield determined by counting the number of events that pass 
all selection criteria and subtracting the background and K® estimates from it, as described in 



B decay mode 


Yield 




43 ± 13 


J/ip7r + ir~ (non-resonant) 


47 ± 15 


-» TT+O 


29 ±6 


J/^vr+vr- (fit) 


90 ± 13 


J/vpir + 7r~ (counting) 


91 ± 17 



where the error is statistical. The results are summarised in Table |T|. 

5 Systematics studies 

The systematic errors on the final branching fraction measurement arise from uncertainties in the 
signal efficiency, fitted yield, number of BB pairs produced and J/ip — > l + t~ branching fraction. 
N B g is known to 1.1% with the dominant contribution to the uncertainty coming from the error 
on the efficiency of the -B-counting selection. B(J/ip — > £ + £~) is known to 1.7% Q. 

The uncertainty on the pion particle identification efficiency is 2.5% per pion. This is determined 
in part by the limited number of events in the control sample. Other contributions are from applying 
different techniques to deduce the kaon contamination in the control sample and deviations between 
the efficiency for pions in the control sample and in the J/ip7r + ir~ sample. Thus the final systematic 
error from pion identification is 5%. Uncertainties on electron and muon particle identification 
efficiencies come from studies of yields in B — > J/ipX data events and give an error of 1.3%. The 
total systematic error on particle identification efficiency is then 5.2%. 

The uncertainty in the determination of the tracking efficiency is 1.3% per track, which sums 
coherently for the four tracks coming from the B°. 

The efficiency of the convergence requirement on the 7r + 7r~ vertex fit has been studied with a 
sample of if>(2S) — > £ + l~ decays. Data and simulation are found to be in good agreement with an 
associated systematic error of 1%. 

The J/ifi p component of the final sample will introduce an uncertainty on efficiency due to the 
unknown polarisation of the p. Studies of efficiency variations between samples of simulated events 
with longitudinal and transverse polarisations lead to a systematic error estimate of 2.5%. 

To account for imperfect simulation of selection variables, the requirements on J/ip invariant 
mass, X^n-K an d are varied within reasonable limits, and the analysis is repeated for each 

selection. The deviation of the efficiency-corrected yield from the fit, in each case, is assigned as a 
systematic error. Varying the size of the J/ip mass interval by ±14 MeV/ c 2 and the | AE\ requirement 
by ±1cj, where a = 13 MeV is the standard deviation observed in data, shows variation at the level 
of 0.4% and 6.8%, respectively. The %3, ^ requirement is varied between 3 and 9 yielding an error 
of 2.1%. These errors are added in quadrature to give the total systematic uncertainty due to 
modelling of the selection requirements. 

To assess the effect of the chosen signal and background shapes on the fitted yields, the fixed 
parameters of these PDFs are varied within their statistical errors, allowing for correlations. This 
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Table 2: Summary of the systematic errors for B(B° — > J/ipir + -K ). 



Source of Uncertainty 


Systematic Error 


B-Counting 


1.1% 


cr(B( J/0 -> ^~)) 


1.7% 


PID efficiencies 


5.2% 


Track efficiencies 


5.2% 


■kit vertex efficiency 


1.0% 


p° polarisation 


2.5% 


Selection variation 


7.1% 


PDF parameter variation 


4.9% 


Simulation statistics 


0.7% 


Total 


11.8% 



produces a total systematic error due to fit parameter variation of 4.9%. The total systematic 
uncertainty from all sources is found to be 11.8%, as summarized in Table [2|. 

6 Cross checks 

The fit was repeated using a simple Breit-Wigner (non-relativistic, width independent of mass) for 
the p° lineshape. While the J/i/j p° yield decreased by 23% and that for J/^7r + 7r~ (non-resonant) 
increased by 27%, the total J/ip7r + Tr~ yield increased by less than 1%. This check is also equivalent 
to varying the parametrization for J/07r + 7r~ (non-resonant) and is strictly a test of the stability of 
the total J/tp 7T + ir~ yield. It should be noted that the simple Breit-Wigner is technically the wrong 
shape to use for the p° parametrization. 

Another way to model the backgrounds is to use a smoothing algorithm on the data rather than 
impose definite PDF shapes. The chosen method creates a Gaussian kernel for each event in order 



to build a shape that describes the input data [11]. The resulting PDF follows fluctuations in the 



input data as accurately as possible and checks how strongly the fitted signal yield depends on the 
chosen method of background parameterization. Changing the background modelling in this way 
alters the total fitted yield by less than one event. 

The signal yield can also be calculated by counting the number of data events passing all the 
cuts and subtracting the estimated numbers of background and J/ip events. This method gives 
91.0 ± 16.6 J/ipir + ir~ events, to be compared to the 89.9 ± 13.3 given by the fit to the M(tt + tc~) 
distribution. 

There are 43.3 ± 12.7 events in the p° resonance component. On the assumption that the non- 
p° signal component is well modelled by J/tpir + 7r~ phase space, we measure B(B° — > J/ifi p°) = 
(2.4 ± 0.7) x 10^ 5 , where the error is statistical only. 



7 Conclusion 

In summary, we have measured the branching fraction for B° meson decay to the final state 
J/ip7r + ir~ to be (5.0 ± 0.7(stat) ± 0.6 (syst)) x 10 -5 . The technique of fitting the mass spectrum 
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of the 7T7T system will become important for future measurements in this channel, particularly as it 
will provide separation between the resonant and non-resonant components given a larger data set. 
It is important to reliably isolate the J/i/j p° component in order to contribute to a measurement of 
the Unitarity Triangle angle 0. 
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